Filler Dependencies of Electrical Conductivity in Nanotube and Nanofiber Composites 
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We report on a model of polymer nanocomposites with fibrous fillers which explicitly considers 
the microscopic filler features and replicates the composites as random distributions of particles 
interconnected via electron tunneling. By exploiting the critical path method, we are able to ob- 
tain simple formulas, applicable to most nanotube and nanofiber composites, which allow to infer 
the overall composite conductivity starting from few parameters like filler volume fraction, size, 
and aspect-ratio. The validity of our formulation is assessed by reinterpreting existing experimen- 
tal results and by extracting the characteristic tunneling length, which is mostly found within its 
expected value range. These results can be used practically to tailor the electrical properties of 
nanocomposites . 
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The inclusion of nanometric conductive fillers such as 
carbon nanotubes [H or nanofibers 0] into insulating 
polymer matrices allows to obtain electrically conductive 
nanocomposites vi^ith unique properties which are widely 
investigated and have several technological applications 
ranging from antistatic coatings to printable electronics 
. A central issue in this domain is to create composites 
with an overall conductivity a which can be controlled 
through the volume fraction 0, the shape of the conduct- 
ing fillers, their dispersion in the polymer matrix, and 
the local inter-particle electrical connectedness. Under- 
standing how these local properties affect the composite 
conductivity is therefore the ultimate goal of any theo- 
retical investigation of such composites. 

A common feature of most random insulator-conductor 
mixtures is the sharp increase of cr once a critical volume 
fraction (j)c of the conductive phase is reached. This tran- 
sition is generally interpreted in the framework of perco- 
lation theory 0,0,01 and associated with the formation 
of a cluster of connected filler particles that spans the 
entire sample. The further increase of a for (f> > (j)c is 
likewise understood as the growing of such a cluster. 

In the case of composites with a polymeric matrix, the 
conducting particles are separated from each other by a 
thin polymeric layer, and the conduction between the 
fillers is originating from quantum mechanical tunnel- 
ing processes 0, H, S ■ Tunneling has consequently 
been implemented in the percolationpicture to explain 
peculiar transport properties near (/)c[E 11, 13, l3- Nev- 
ertheless, although a good understanding of the factors 
affecting (j),, exists, a formulation of cr((^) for (j) > <j)c with 
predictive power is still missing. This especially applies 
to polymer composites with fillers such as nanotubes and 
nanofibers which have diameters which are comparable 
with the typical distances of tunneling. In these cases, 
the detailed behavior of tunneling is expected to take an 
explicit role in a{(t)) and should be taken into account. 

In this Letter we address this problem by assuming 
that the conducting fillers form a network of globally 



connected sites via tunneling processes [31 ■ We show 
that the calculated conductivity of composites with dif- 
ferent filler aspect-ratios can be re-obtained by the criti- 
cal path method, allowing for a simple analytical formula 
for cr{(j}), valid for many regimes of interest, and having 
an explicit dependence on microscopic quantities such 
as the geometrical dimensions of the fillers and the tun- 
neling length. We further show that, by reinterpreting 
conductivity data of real nanotubes and nanofibers com- 
posites, our formulation leads to consistent estimations 
of the tunneling length for these systems. 

In modeling the conductor-insulator composite mor- 
phology, we treat the conducting fillers as identical im- 
penetrable objects dispersed in a continuous insulating 
medium, with no interactions between the conducting 
and insulating phases. Furthermore, in order to describe 
rod-like filler particle shapes, we consider impenetrable 
spheroids (ellipsoids of revolution) up to the the extreme 
prolate [a/b ^ 1) limit, where a and h are the spheroid 
polar and equatorial semi-axes respectively. Distribu- 
tions of spheroids with different aspect-ratios a/b and 
volume fractions </> were generated inside a cubic cell with 
periodic boundary conditions via random sequential ad- 
dition [TBI and successively relaxed through monte-carlo 
(MC) sweeps [l3,[lB|- High density configurations were 
obtained by combining MC sweeps with particle inflation 
[l6| . The absence of global oricntational order was 
verified using the nematic order parameter [l7j . 

In considering the overall conductivity arising in such 
composites, we assume the spheroids of our distributions 
to be perfectly conductive and attribute to each spheroid 
pair i, J a tunneling conductance of the form Q 



9ij = go exp 



25^ 



(1) 



where ^ is the characteristic tunneling length, which mea- 
sures the electron wave function decay within the poly- 
mer, and Sij is the minimal distance between the two 
spheroid surfaces. For spheres (a/6 = 1), Sij is simply the 
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FIG. 1: (Color online) (a) Volume fraction </!> dependence of the tunneling conductivity a for a system of aspect-ratio a/h = 10 
hard prolate spheroids with different characteristic tunneling distances ^/D with D — 2a. Results from Eq. ([2| (with cro = 0.179) 
are displayed by dotted lines, (b) Tunneling conductivity in a system of hard prolate spheroids with different aspect-ratios a/b 
for the same £,/ D. Dotted lines; results from Eq. ([2]) with (Tq = 0.124 for a/h = 2, and ctq = 0.115 for a/h = 1. Inset: schematic 
illustration of the tunneling conductivity crossover for the cases a/b — 1, 2, and 10. 



center-to-center distance minus twice the radius of the 
two spheres, while for the general case {a/b ^ 1). 5ij de- 
pends also on the relative orientation of the spheroids and 
can be obtained from a numerical procedure described in 
Ref. [3| • The specific value of the pre-factor go may vary 
for different composites and, compared to the exponential 
term in Eq. ([1]), has a weaker dependence on Sij, which 
we neglect by treating go as a constant with unit value. 
We also note that in writing Eg. ([T]) we disregarded par- 
ticle charging effects [3 which, at constant room 
temperature, that is the case considered here, and if the 
particles arc identical and sufficiently large, may be ne- 
glected compared to the tunneling decay. 

To evaluate the global conductivity of the system, the 
full set of bond conductances given by equation Eq. ^ is 
mapped on a resistor network and the overall conductiv- 
ity is calculated through numerical decimation of the re- 



sistor network [13|, |20| . To reduce computational times to 
manageable limits, a cutoff distance Sco is introduced in 
order to reject negligibly small bond conductances. This 
was typically equal to 4a for low density configurations, 
while, for spheroids at high densities, it was reduced to 
Sco = 2a for a/b = 1, 2 and Sco = a for a/6 = 10. 

In Fig. [ija) we report the so-obtained conductivity a 
values (symbols) as function of the volume fraction of 
prolate spheroids with aspect-ratio a/b ~ 10, and for 
different values of £,/D, where D = 2a. Each symbol 
is the outcome of Nn ~ 200 realizations of a system 
of Np ^ 1000 spheroids. The logarithm average of the 
results was considered since, due to the exponential de- 
pendence of ([T]); the distribution of the computed con- 
ductivities was approximately of log-normal form. The 



strong reduction of a for decreasing is a direct conse- 
quence of the fact that as (f> is reduced, the inter-particle 
distances get larger, leading in turn to a reduction of 
the local tunneling conductances ([T]). Furthermore, as 
shown in Fig. W[h), such reduction depends strongly on 
the aspect-ratio of the conducting fillers: as a/b increases, 
the composite conductivity drops for much lower values 
of <f) for fixed ^. However, in real composites, we have 
to take into account the polymer matrix whose intrin- 
sic conductivity (t„i, which falls typically in the range 
(Tm ~ lO^^'^ ^ 10"^* S/cm, prevents an indefinite drop 
of a. This is schematically illustrated in the inset of 
Fig. [Hb). Now, if we identify (pc with the volume frac- 
tion at which a ~ a^n, then fillers with larger aspect- 
ratios entail lower values of 0c, consistently with what is 
commonly observed. This also implies that a model of 
the composite conductivity based on the tunneling con- 
tribution of Eq. IT]) alone will be representative only if (jj 
is sufficiently larger than 4>c to consider the effect of am 
negfigible. 



We are now going to show that the strong dependence 



of cr(0) on a/b and ^ of Fig. 
critical path (CP) method 



2] can be reproduced by the 
2ll . [2^ applied to our sys- 



tem of impenetrable spheroids. For the tunneling con- 
ductances of Eq. ID), this method amounts to keep only 
the subset of conductances gij having Sij < Sc, where 
Sc, which in turn defines the characteristic conductance 
9c = <?o exp(— 2(5c/'C), is the largest among the Sij dis- 
tances such that the so-defined subnetwork forms a con- 
ducting cluster spanning the sample. Next, by assigning 
gc to all the (larger) conductances of the subnetwork, a 
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FIG. 2: (Color online) Re-scaled critical distances 5c/ D versus 
for the prolate spheroids with a/b — 1, 2, 10, 20, and 100, 
and for the impenetrable spherocylinders of Refs. [l7l.[2^. For 
a/b =1 our results are plotted together with those of Ref. [2^ . 
The solid line is the asymptotic behavior for a/6 ^ 1. 



CP approximation for a is 



(T ~ ctq exp 



2(5,(0, Q, 6) 



(2) 



where ctq is a pre-factor proportional to go- The signif- 
icance of Eq. ([2]) is that it reduces the conductivity of 
a distribution of hard objects electrically connected by 
tunneling to the computation of the geometrical "criti- 
cal" distance 5c- In practice, 5c can be obtained by coat- 
ing each impenetrable spheroid of a distribution with a 
penetrable shell of constant thickness (5/2, and by consid- 
ering two spheroids as connected if their shells overlap. 
5c is then the minimum value of 5 such that, for a given 
0, a cluster of connected spheroids spans the sample. In 
order to illustrate the efficacy of the CP method, we have 
calculated 5c for the same parameters of Fig. [1] (see be- 
low), obtaining an excellent agreement of Eq. ([2]) (dotted 
lines) with the outcomes of the full numerical decimation 
of the resistor network. 

We shall now illustrate that for sufficiently elongated 
prolate spheroids, a simple relation exists that allows to 
estimate 5c with good accuracy. In virtue of Eq. ([2]) this 
means that we can formulate explicit relations between 
a and the shapes and concentration of the conducting 
ffllers. 

To extract 5c we followed the route outlined in 
Ref. [l3|- For given combinations of spheroid aspect- 
ratios and distance 5c, we searched the critical (percola- 
tion) volume fraction by recording for several (j) the prob- 
ability of having a cluster spanning the opposite ends of 
the simulation cell, and identified it with the concentra- 
tion for which this probability is equal to 1/2. For every 



configuration, we used = 40 for the smallest values 
of 5c up to Nn = 500 for the largest ones. The parti- 
cle number varied between Np ^ 2000 (a/6 = 1,2) and 
Np ^ 8000 (a/5 = 100). Relative errors on the per- 
colation volume fraction were in the range of a few per 
thousand. 

In Fig. [2] we report the calculated values of 5c/ D as a 
function of volume fraction (f) for spheres (a/6 = 1, to- 
gether with the results of Ref. [31), and for a/b = 2, 
10, 20, and 100. In the figure, 5c/ D is multiplied by the 

ratio Vsphere/Vobject = (ci/b)^, where Vgphere = ttD^/6 

is the volume of a sphere with diameter equal to the 
major axis of the prolate spheroid and Vobject is the vol- 
ume of the spheroid itself. For comparison, we plot in 
Fig. [2] also the results for impenetrable spherocylinders 
of Rcfs. [TtI . 24 1 . These are formed by cylinders of ra- 



dius R and length L, capped by hemispheres of radius R, 
so that a — R + L/2 and b = R, and Vgphere /Vobject = 
(a/6)V[(3/2)(a/6) - 2] ~ {2/3){a/b)^ for a/b > 1. As it 
is apparent, for sufficiently large values of a/6 the simple 
re-scaling transformation collapses both spheroids and 
spherocylinders data into a single curve. This holds true 
as long as the aspect-ratio of the spheroid plus the pen- 
etrable shell (a -I- (5c/2) /{b + 5c/2) is larger than about 5. 
In addition, for (p < 0.03 the collapsed data are well ap- 
proximated by 5c{Vsphere /Vobject) / D = OA/ (j) (solid line 
in Fig. (2), leading to the following asymptotic formula: 



5c/ D 



lib/af 



(3) 



where 7 = 0.4 for spheroids and 7 = 0.6 for spherocylin- 
ders. Equation ([3]) is fully consistent with the scaling law 
of Ref. [25| based on the second- virial approximation for 
semi-penetrable spherocylinders. 

We are now in the position to evaluate the tunneling 
driven conductivity of random distributions of prolate 
objects for a > <Tm- By substituting Eq. ^ into Eq. ([2|) 
we obtain 



(T ~ (To exp 



2D7(6/a) 
"1 ~ 



(4) 



From the previously discussed conditions on the validity 
of the asymptotic formulas for 5c/ D it follows that the 
above equation will hold when {b/aY ^ </> 0.03. If 
these conditions are not met, values of 5c can be anyway 
extrapolated from the data of Figs. [2] 

Let us now show how the above outlined formalism 
may be used to re-interpret the experimental data on 
the conductivity of different nanocomposites with fibrous 
ffilers that can be found in the literature. In Fig. [3ja) we 
report measured data of ln((T) versus (j) for polymer com- 
posites filled with Cu nanofibers [2^, and carbon nan- 
otubes [131 • Equation ^ implies that the same data can 
be more profitably re-plotted as a function of 5c, instead 
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FIG. 3: (Color online) (a) Natural logarithm of the conductiv- 
ity (J as a function of the volume fraction (ji for Cu nanofibers- 
polystyrene[2^, and single-wall carbon nanotubes-epoxv[27| 
composites. When, for a given concentration, more then one 
value of a was given, the average of In a was considered, (b) 
The same data of (a) re-plotted as function of the correspond- 
ing critical distance Sc. Solid lines are fits to Eq. ([5}. (c) 
Characteristic tunneling distance ^ values for different poly- 
mer nanocomposites as extracted by means of Eq. ((5]) 



of d). Indeed, from 



In(CT) = --(5c + In(o-o), 



(5) 



a linear behavior with slope —2/^ is expected, indepen- 
dently of the specific value of ctq, which allows for a direct 
evaluation of the characteristic tunneling distance ^. The 
exact value of ^ is a much investigated topic [il], and is 
expected to be between a fraction of nm and a few nm 
d, m, m, [lOl. By using the values of D and a/b 
provided in Refs. (2^, H^] and our formula ^ for Sc, 
we find indeed an approximately linear dependence on 
[Fig. [3jb)], from which we extract 1.65 nm for the 
nanotubes, and 5.9 nm for the nanofibers. 

We further applied this procedure to several published 
data on polymer-based nanotube [H and nanofiber Q 
composites with fillers having a/b ranging from ~ 20 up 
to ~ 2 • 10^. The results are collected in Fig.[3Jc), show- 
ing that most of the so-obtained values of the tunneling 
length ^ are comprised between ~ 0.1 nm and ~ 10 nm. 
This is a striking result considering the number of fac- 
tors that make a real composite deviate from an idealized 
model. Most notably, fillers may have non-uniform size, 
aspect-ratio, and geometry, they may be oriented, bent 
and/or coiled, and interactions with the polymer may 
lead to agglomeration, segregation, and sedimentation. 



Furthermore, composite processing can alter the prop- 
erties of the pristine fillers, e.g. nanotube or nanofiber 
breaking [which may explain the downward drift of ^ for 
high aspect-ratios in Fig. [3I^c)]. In principle, deviations 
from ideality can be included in the present formalism by 
evaluating their effect on 5c {2^. It is however interest- 
ing to notice that all these factors have often competing 
effects in raising or lowering the composite conductivity, 
and Fig. [3^c) suggests that in average they compensate 
each other to some extent, allowing tunneling conduction 
to emerge strongly from a{(j)) as a visible characteristic 
of nanocomposites. 

In summary, by extending the critical path method, 
we have mapped the tunneling driven conductivity cr of 
dispersions of fibrous conducting fillers into a geometri- 
cal percolation problem, which has permitted us to for- 
mulate the filler dependencies of a in terms of a single 
parameter: the critical distance Sc- We have shown that 
for sufficiently large aspect-ratios of the fillers, Sc, and 
so cr, can be expressed by a simple formula applicable to 
real dispersions of nanotubes and nanofibers. To validate 
our formulation, we have analyzed published conductiv- 
ity data for several nanofibers and nanotube composites 
and extracted the corresponding values of the tunneling 
length ^, which is mostly found within its expected range. 
The above outlined procedures can be likewise used as 
guidelines to tailor the electrical properties of real com- 
posites, and can be generalized to include different filler 
shapes, filler size and/or aspect-ratio polydispersity, and 
interactions with the polymer. 
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